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We present a search for point sources of high energy neutrinos using 3.8 years of data recorded by AMANDA- 
II during 2000-2006. After reconstructing muon tracks and applying selection criteria designed to optimally 
retain neutrino-induced events originating in the Northern Sky, we arrive at a sample of 6595 candidate events, 
predominantly from atmospheric neutrinos with primary energy 100 GeV to 8 TeV. Our search of this sample 
reveals no indications of a neutrino point source. We place the most stringent limits to date on E~ 2 neutrino 
fluxes from points in the Northern Sky, with an average upper limit of E 2 & v + Vt < 5.2 x 10 -11 TeV cm 
s _1 on the sum of and v T fluxes, assumed equal, over the energy range from 1.9 TeV to 2.5 PeV. 



PACS numbers: 95.85.Ry, 95.55.Vj, 98.70.Sa 

I. INTRODUCTION 

Detecting extraterrestrial sources of high energy (>TeV) 
neutrinos is a longstanding goal of astrophysics. Neutrinos 
are neither deflected by magnetic fields nor significantly at- 
tenuated by matter and radiation en route to Earth, thus neu- 
trino astronomy offers an undistorted view deep into the high 
energy universe. Particularly, neutrinos offer an opportunity 
to probe the sources of high energy cosmic rays, which re- 
main unknown. Potential cosmic ray sources include galactic 
microquasars and supernova remnants as well as extragalactic 
sources such as active galactic nuclei and gamma ray bursts. 
These objects are thought to accelerate protons and nuclei in 
shock fronts via the Fermi mechanism yl, resulting in power 
law energy spectra E a , with a ~ —2. A fraction of the 
energized particles interact with local matter and radiation, 
producing pions. The neutral pions decay into high energy 
photons, and the charged pions ultimately produce neutrinos 
with a flavor ratio v e :v^:v T ^1:2:0, mixing to approximately 
1:1:1 at Earth because of vacuum flavor oscillations. Obser- 
vations of TeV gamma rays [2J, yl |4J] hint at possible cosmic 
ray source locations but currently cannot separate neutral pion 
decay spectra from inverse Compton emission. The Auger 
collaboration has reported a correlation of arrival directions 
of the highest energy cosmic rays with active galactic nu- 
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clei |Ht]; however, a similar correlation has not been observed 
by HiRes ((J. Identification of a high energy neutrino point 
source would provide an unambiguous signature of energetic 
hadrons and cosmic ray acceleration. Neutrino flux predic- 
tions exist for many potential sources flSEEEllllIlIlIllj but 
no high energy neutrino point source has yet been identified 
IB HQ. 

The search for high energy neutrino point sources is a ma- 
jor objective of the Antarctic Muon And Neutrino Detector 
Array (AMANDA). High energy leptons are produced in the 
Earth by charged-current neutrino interactions. In transparent 
matter, a cone of Cherenkov photons propagates from the lep- 
ton track according to the optical properties of the medium. 
AMANDA-II is an optical Cherenkov detector consisting of 
677 optical modules arranged in 19 strings frozen ~ 1500 m to 
^2000 m deep in the ice sheet at the geographic South Pole. 
Approximately 540 modules in the core of the array showing 
stable performance are used in this search. Each module con- 
tains a 20 cm diameter photomultiplier tube (PMT) optically 
coupled to an outer glass high-pressure sphere. PMT pulses 
are propagated to surface electronics, and, when the trigger 
threshold of 24 discriminator crossings ("hits") within 2.5 /is 
is satisfied, the pulse leading edge times are recorded. The 
leading edge times along with known detector geometry and 
optical properties of South Pole ice [16] allow reconstruction 
of tracks passing through the detector 11711 . High energy elec- 
trons produce short electromagnetic cascades with little direc- 
tional information of the primary neutrino. Muons produced 
in the ice and bedrock, on the other hand, propagate up to 
several kilometers to the detector and their tracks are recon- 
structed with 1.5°-2.5° median accuracy depending on en- 
ergy and zenith angle. Tau leptons decay rapidly and produce 
tracks too short for reconstruction below ^PeV energies. Tau 
decay, however, contributes high energy muons with a branch- 
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Year Livetime Total Events Filtered Events Final Selection 



2000 


197 d 


1.37 xlO M 


1.63 xl0 b 


596 


onn i 
2UU1 


1 nQ ^ 

i93 a 


o nn w i n9 


1 nn w i n6 
i.9UX 1U 


854 


2002 


204 d 


1.91 xlO 9 


2.10xl0 6 


1009 


2003 


213 d 


1.86xl0 9 


2.22x10 s 


1069 


2004 


194 d 


1.72xl0 9 


2.09 xlO 6 


998 


2005 


199 d 


2.06 xlO 9 


5.21 xlO 6 


1019 


2006 


187 d 


2.00 xlO 9 


4.89 xlO 6 


1050 


Total 


1387 d 


12.92xl0 9 


20.04 xlO 6 


6595 



TABLE I: AMANDA livetime and event totals. 



ing ratio of 17.7% lfl3"l[l8ll . and these muon tracks can be re- 
constructed. We thus search for upward propagating muons 
produced in the Earth by (D^) and v T (D T ) fluxes following 
roughly an E~ 2 energy spectrum. While downward neutrino 
induced muons also trigger the detector, such events are diffi- 
cult to distinguish from downward muons produced by cosmic 
ray air showers. Located at the South Pole, AMANDA-II is 
thus most sensitive to neutrino fluxes from the Northern Sky. 
Air showers also produce neutrinos, and this atmospheric neu- 
trino flux lfl9il20ll is the main background for our search. 

Here we present the results of a search for astrophysical 
point sources of high energy neutrinos using 3.8 years of data 
recorded by AMANDA-II during 2000-2006, extending the 
previous five-year analysis lfl3ll with data from the final two 
years of standalone operation and improving our sensitivity 
by a factor of ~2. We report flux limits for a catalog of 26 
selected source candidates along with results of a search for 
neutrino sources over the entire Northern Sky. Additionally, 
we report results from a search for neutrino emission from 
gamma ray sources identified by Milagro [2] and a search for 
event angular correlations. In all cases, we observe no indica- 
tions of an astrophysical neutrino point source. 

II. DATA SELECTION 
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FIG. 1 : Zenith angle (9) distributions for data and simulation at sev- 
eral reduction levels. Reconstructed (solid) and true (fine dotted) 
zenith angle distributions are shown for CORSIKA 1 23] cosmic ray 
muon simulation at retrigger level, and reconstructed zenith angle 
distributions are shown for atmospheric neutrino simulation (dotted) 
and data (circles) at retrigger level, filter level, and final selection. 
We also show the reconstructed zenith angle distribution of a diffuse 
E~ 2 neutrino flux at the current limit 121] using our final selection 
(dash-dotted). 



As illustrated in Fig. Q] AMANDA-II records O(10 9 ) 
events per year from downward propagating muons produced 
by cosmic ray air showers, O(10 3 ) events per year from at- 
mospheric neutrinos, and O(10) high quality events per year 
from astrophysical E~ 2 neutrino fluxes given current limits 
j2lll . We attempt to isolate these neutrino events from the 
downward muon background in a computationally efficient 
manner. We exclude data taken during periods of detector in- 
stability and significant maintenance, which include the aus- 
tral summer (November 1 through February 15). After ac- 
counting for deadtime in data acquisition electronics, nomi- 
nally ~15% of uptime, we have accumulated 1387 days (3.8 
years) of livetime with 1.29xl0 10 events during seven years 
of operation (table [I]). 

Events are first processed to remove hits induced by electri- 
cal cross talk, hits from unstable modules, and isolated noise 
hits fl7tl . and events which no longer pass the trigger criteria 
are discarded. These retriggered events are then reconstructed 
with the fast pattern matching algorithms DirectWalk (DW) 



lfl7ll and JAMS [22] which identify muon tracks within events. 
For optimal efficiency, our upgoing event selection requires 
both zenith angles 9dw and 9 jams greater than 70°-80°. 

Two CPU intensive maximum likelihood reconstructions 
are applied to events passing the above selection criteria. First, 
we apply an unbiased likelihood (UL) fit seeded with the Di- 
rectWalk and JAMS reconstructed tracks and 30 additional 
random track directions. The likelihood function [ 17] param- 
eterizes the probability of observing the obtained geometry 
and leading edge times of hit modules in terms of track zenith 
angle, azimuthal angle, and position. The likelihood is max- 
imized with respect to these parameters (in practice, the neg- 
ative logarithm of the likelihood is numerically minimized), 
yielding the best fit track zenith and azimuthal angles, and the 
fit result from the seed yielding the maximum likelihood is 
chosen as the reconstructed track. A 64 seed Bayesian likeli- 
hood (BL) fit is also done, using the downgoing muon zenith 
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FIG. 2: Distributions of data and atmospheric neutrinos at filter level and final selection level for several parameters and zenith angles 8 > 95° 
(top and left), and zenith angle distribution for the selected 6595 neutrino candidate events compared with model predictions fl9l [20ll for 
atmospheric neutrinos (bottom right). 



angle distribution as a Bayesian prior. With the additional cut 
&ul > 80°, our upgoing event filter reduces the downward 
muon background by a factor of ^650 relative to trigger level 
(table©. 

After this cut, O(10 6 ) misreconstructed downward muon 
events per year remain, which still outnumber atmospheric 
neutrinos by roughly three orders of magnitude. The vast ma- 
jority of these events are removed by the following four topo- 
logical parameters, shown in Fig. [2] 

• The likelihood ratio of the UL and BL fits. Downgoing 
muon background events misreconstructed as upgoing 
by the UL fit typically are also fit well with the down- 
ward biased BL reconstruction, whereas true upgoing 
events are not. Therefore, the UL/BL likelihood ratio 



tends to be higher for upgoing events. 

• The angular uncertainty of the UL fit, described further 
in section III. Misreconstructed events generally have 
large angular uncertainty. 

• The smoothness, or homogeneity of the hit distribution 
along the UL track [17]. High quality events contain 
photon hits along the entire length of the track and have 
smoothness values near zero, whereas hits from misre- 
constructed events tend to distribute near the beginning 
or end of the track and have smoothness values near +1 
and —1, respectively. 

• The UL track direct length, obtained by projecting di- 
rect hits backward to the UL track at the Cherenkov 
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5=90° 




FIG. 3: Equatorial sky map of 6595 events recorded by AMANDA-II from 2000-2006. A table of the events is available [26 



angle and taking the distance along the track between 
the first and last. We select direct hits, compatible 
with relatively unscattered photons and arriving on- 
time with the Cherenkov cone, using the time window 
— 15 ns < t — t c h < 25 ns 11711 . Hits from misrecon- 
structed events rarely follow the muon-Cherenkov tim- 
ing pattern over significant distances, resulting in short 
lengths. 

For the zenith angle region 91.5° < 8 < 180° we use the fol- 
lowing zenith angle dependent cuts, optimized to yield maxi- 
mum sensitivity 02411 : 

\og{UL/BL) > 34 - 25 • $(cos 6 + 0.15) 

Oi <C 3.2 — 4 • $(- cos 9 - 0.75) 
|Smoothness| < 0.36. 

Here = x for positive x, and $(x) = for x < 0. We 
use a support vector machine (SVM) 112511 trained on the four 
parameters to improve event selection in the near-horizontal 
region 80° < 8 < 91.5°. Events with SVM quality of zero or 
less are consistent with misreconstructed muon background, 
while events with larger values of SVM quality are increas- 
ingly consistent with quality muons. We apply the cut: 

SVM Quality > 1 - 12 • $(cos 8 - 0.023) 

Application of these quality cuts yields 6595 neutrino candi- 
date events [26] (Fig. [3j. 

Simulations of two atmospheric neutrino flux models lfl9[ 
Eoll . with events generated by ANIS l27ll an d resultant muons 
propagated to the detector with MMC [28], both agree with 
data in track quality parameter distributions and zenith angle 
(Fig. |2]i within the ^30% uncertainty in these flux predic- 
tions. Application of the filter selection and final quality cuts 
to this simulation yields an atmospheric neutrino efficiency 
of 30% relative to retrigger level for 8 > 90°. The contri- 
bution of misreconstructed downward muons has been esti- 
mated by subtracting the simulated atmospheric neutrino rate, 
after renormalizing it for a more stringent selection yielding 
a nearly pure neutrino sample. The muon contamination has 
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FIG. 4: Effective area for averaged and (solid) and averaged 
v T and v-r (dashed) neutrino fluxes for several declination ranges. 



been found to be less than 5% for 8 > 95° (declination 5 > 
5°), but the contamination is more significant near the equator 
and dominates events in the Southern Sky. A parallel analysis 
of these atmospheric neutrino events has revealed no evidence 
of new physics such as violation of Lorentz invariance and 
quantum decoherence M29I1 . We simulate and v T events 
from 10 GeV to 100 PeV with an identical software chain, 
and this simulation is used to calculate neutrino effective area, 
shown in Fig. |4j and flux limits for neutrino sources with 
E~ 2 energy spectra. The central 90% of such signal events 
fall within the energy range 1.9 TeV to 2.5 PeV. The median 
accuracy of the UL fit when applied to simulated events fol- 
lowing an E~ 2 energy spectrum is 1 .5°-2.5°, shown in Fig. [5] 
The absolute pointing accuracy of AMANDA has been con- 
firmed by observing downgoing muon events coincident with 
well-reconstructed air showers recorded by SPASE IU7I1 and 
events coincident with IceCube. 
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FIG. 5: Angular deviation between neutrino and UL fit track for sim- 
ulated E~ 2 muon neutrino events from several declination ranges. 



FIG. 7: Distributions of angular deviation between true and recon- 
structed tracks for simulated neutrino events over several ranges of 
estimated angular uncertainty. 



"[Atmospheric Neutrino Flux 




Log io (E v /GeV) 

FIG. 6: Energy distribution of events passing selection criteria for 
simulated atmospheric neutrino background 1 19j in a 3.5° bin and an 
E~' 2 point source with flux § Vfl+VT = 10" 10 TeV cm~ 2 s _1 . Such 
a source would be detected at 5<r in approximately 40% of trials. 



III. SEARCH METHOD 



The remaining background, mostly atmospheric neutrinos, 
is difficult to reduce further without significantly decreasing 
signal efficiency. Neutrinos from E~ 2 sources are typically 
more energetic than atmospheric neutrinos (Fig. [6]), which fol- 
low a steeper ^E~ 3 7 energy spectrum. We search our sam- 
ple of 6595 events for excesses above the atmospheric neu- 
trino background both in direction and event energy using an 
unbinned maximum likelihood search method [30], providing 
direction and energy discrimination on an event-by-event ba- 
sis by incorporating an event angular resolution estimate and 
energy estimate. 



A. Event Angular Uncertainty Estimation 

Our ability to reconstruct muon tracks in AMANDA par- 
tially depends on event topology. A muon track passing 
through a larger portion of the detector or giving hits in a 
larger number of modules should, on average, reconstruct 
with better angular resolution due to a longer lever-arm or 
larger number of measurement points, respectively. We there- 
fore estimate the resolution of each UL track by evaluating 
the likelihood space near the maximum 113 ill . As the track 
zenith angle and azimuthal angle coordinates (8, </>) move 
away from the best fit track values (6, </>), the quantity log C 
decreases parabolically from its maximum. The likelihood 
ratio —2 • log ( C [ S ~'^l ) is evaluated on a grid of zenith and az- 
imuthal angles near the best track, and the resulting values are 
fit to a paraboloid with the form 




where the x and y axes are fit and do not necessarily cor- 
respond to zenith and azimuthal angles. The two errors a x 
and <j y are then geometrically averaged into a single, circu- 
lar error crj. The paraboloid fit is thus a convenient approxi- 
mation of the likelihood space, reducing the complex map of 
—2 • log into just Cj, The corresponding spatial prob- 

ability density estimate at an angular distance ^ is then: 




Distributions of the angular deviation between true and recon- 
structed neutrino tracks for several ranges of estimated an- 
gular uncertainty (Fig. |7J show the correlation between es- 
timated angular uncertainty and track reconstruction error. 
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FIG. 8: Muon energy distributions for four ranges of N c h (top), and 
simulated N c u distributions for atmospheric neutrinos 1 19] and E -2 , 
E -2 ' 5 , and E -3 power law neutrino spectra (bottom). 



B. Event Energy Estimation 



The amount of light deposited in the detector depends 
strongly on muon energy above ^1 TeV, and thus the num- 
ber of hit modules (N c h) provides an approximate measure of 
event energy. Distributions of muon energy for several ranges 
of N c h (Fig. [8j show the performance of N c h as a muon en- 
ergy estimator, with a ler uncertainty in \og 10 (E^/GeV) of 
0.65. Rather than measure event absolute energy, it is more 
relevant for a neutrino search to assess the compatibility of an 
event with expected astrophysical neutrino spectra, assumed 
to follow a power law. From simulations, we tabulate N c h 
probabilities for spectral indices 1 < 7 < 4 in bins of 0.01 
and for atmospheric neutrinos 1 19], shown in Fig. [8] This ta- 
ble yields the probability of observing a given N c h value from 
a source with a power law energy spectrum relative to observ- 
ing the value from background atmospheric neutrinos. 



C. Maximum Likelihood Method 

For a source with position x s , giving n s events against a 
background of N — n s events, the probability density is 



N 



(3) 



where S and B are the probability densities for signal and 
background, respectively. The likelihood function is 



/V 



N 



-S, 



(4) 



where i runs over the selected events. Events are assumed 
to have an angular error distributed according to a Gaussian 
given by the event angular uncertainty a i7 and signal events 
are assumed to follow a power law energy spectrum with spec- 
tral index 7. The signal probability density for an event at Xi 
is 



2ira? 



(5) 



where \xi — x s \ is the angular distance between the event and 
assumed source position. In practice, we only include events 
with declinations ±8° of the source declination since events 
outside this band have extremely low signal probabilities, and 
we set N to be the number of events in this declination band. 
The background probability over this band is roughly constant 
and given by 



B, 



P(N ch!i \<p atm ) 



(6) 



The likelihood L is maximized (again, — log L is numerically 
minimized) with respect to n s and 7, resulting in best fit signal 
strength h s and spectral index 7. The data are then compared 
to the null, background-only hypothesis (n s — 0) to determine 
relative compatibility. We use as our test statistic 



A 



f£(n s =Q)\ 



(7) 



Larger values of A reject the null hypothesis with increasing 
confidence, shown in Fig. (O. The significance of a particu- 
lar value of A is determined by comparing the obtained value 
to the distribution of test statistic values at the same location 
from data randomized in right ascension, and we denote as p 
the fraction of randomized data sets with higher test statistic 
values. This method, by using unbinned event-by-event en- 
ergy and directional discrimination, improves the sensitivity 
to E~ 2 neutrino fluxes by more than 30% relative to the pre- 
vious method ll32ll using angular bins. 



IV. SEARCH FOR POINT SOURCES IN THE NORTHERN 

SKY 

We first apply the search to a predefined list of 26 ener- 
getic galactic and extragalactic objects, including many TeV 
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FIG. 9: Integral distribution of the test statistic for background at 
5=42.5° with 3a and 5a thresholds indicated and statistical uncer- 
tainty shaded in gray (top), and distribution of the test statistic for 
background and 6, 12, and 18 added E~ 2 signal events at 5=42.5° 
(bottom). 



gamma ray sources. For each source location, we compute the 
value of the unbinned search test statistic A. Flux upper limits 
are computed from the test statistic using Feldman-Cousins 
unified ordering [33]. Systematic uncertainties are incorpo- 
rated into the limit calculation using the method of Conrad 
et al. [34] as modified by Hill ll35ll We estimate the total 
systematic uncertainty in our event rate expectations for E~ 2 
fluxes to be 17%, summarized in table [TTJ Significant con- 
tributions include the absolute sensitivity of optical modules 
(9%), neutrino interaction cross section (8%), bias in event se- 
lection between data and simulation (7%), and photon propa- 
gation in the ice (5%), determined by detailed detector stud- 
ies presented in [13]. Additionally, we evaluate bias in re- 
construction accuracy by comparing distributions of event an- 
gular resolution estimates (<Ti) with those from point source 
simulations. We find the angular resolution estimates in sim- 
ulation are typically 8% smaller, and adjusting our simulated 
point spread by this factor results in flux limits 7% higher. 
Other known sources of systematic uncertainty, including un- 
certainties in optical module timing resolution and the search 
method, total less than 4%. Limits on v u + v T fluxes at 90% 



Neutrino cross section and rock density 


±8 


% 


Optical module sensitivity 


+2 
-9 


% 


Photon propagation 


±5 


% 


Event selection bias 


+0 
-7 


% 


Event reconstruction bias 


+0 


% 


Other known sources 


< 4 


% 


Total 


+10 
-17 


% 



TABLE II: Systematic errors in event rate expectations for point 



sources with E energy spectra. 



Candidate 


<5(°) 


a(h) 


6 Qn 
^ yu 




*(°) N 


3C 273 


2.05 


12.49 


8.71 


0.086 


9 1 
L. 1 




SS 433 


4.98 


19.19 


3.21 


0.64 


1 1 


i 


GRS 1915+105 


10.95 


19.25 


7.76 


0.11 


Z.J 


Q 
O 


M87 


12.39 


12.51 


4.49 


0.43 




~2 
J 


PKS 0528+134 


13.53 


5.52 


3.26 


0.64 




n 
\t 


3C 454 3 


16.15 


22.90 


2.58 


0.73 


Z.J 


c 

J 


fjp m i n rro 
VJClllllliid. 


1 7 77 


(S S7 

U.J / 


12.77 


OOXfi 


2.3 


2 


Crab Nebula 


22.01 


5.58 


9.27 


0.10 


2.3 


1 


GRO J0422+32 


32.91 


4.36 


2.75 


0.76 


2.2 


3 


Cyg X-l 


35.20 


19.97 


4.00 


0.57 


2.1 


3 


MGRO J20 19+37 


36.83 


20.32 


9.67 


0.077 


2.1 


7 


4C 38.41 


38.14 


16.59 


2.20 


0.85 


2.1 


4 


Mrk421 


38.21 


11.07 


2.54 


0.82 


2.1 


3 


Mrk 501 


39.76 


16.90 


7.28 


0.22 


2.0 


6 


Cyg A 


40.73 


19.99 


9.24 


0.095 


2.0 


3 


Cyg X-3 


40.96 


20.54 


6.59 


0.29 


2.0 


8 


Cyg OB2 


41.32 


20.55 


6.39 


0.30 


2.0 


8 


NGC 1275 


41.51 


3.33 


4.50 


0.47 


2.0 


4 


BL Lac 


42.28 


22.05 


5.13 


0.38 


2.0 


2 


H 1426+428 


42.68 


14.48 


5.68 


0.36 


2.0 


3 


3C66A 


43.04 


2.38 


8.06 


0.18 


2.0 


6 


XTE Jill 8+480 


48.04 


11.30 


5.17 


0.50 


1.8 


3 


1ES 2344+514 


51.71 


23.78 


5.74 


0.44 


1.7 


2 


Cas A 


58.82 


23.39 


3.83 


0.67 


1.6 


2 


LS I +61 303 


61.23 


2.68 


14.74 


0.034 


1.5 


5 


1ES 1959+650 


65.15 


20.00 


6.76 


0.44 


1.5 


5 



TABLE III: Flux upper limits for 26 neutrino source candidates: 
Source declination, right ascension, 90% confidence level upper lim- 
its for ^ M + v T fluxes with E~ 2 spectra (_E 2 $„ fl +„ T < $90 x 
10" 11 TeVcm" 2 s -1 ) over the energy range 1.9 TeV to 2.5 PeV, 
pre-trials significance, median angular resolution of primary neu- 
trino, and number of events inside a cone centered on the source 
location with radius equal to the median point spread. Since event 
energy is an important factor in the analysis, the number of nearby 
events does not directly correlate with pre-trials significance. 



confidence level and chance probabilities (p) are shown in ta- 
ble |nil Limits on fluxes alone correspond to half these 
values. The highest significance is found for Geminga with p 
= 0.0086. The probability of obtaining p < 0.0086 by chance 
for at least one of 26 sources is 20% and is therefore not sig- 
nificant. 

We then apply the search to declinations —5° < S < 83° 
on a 0.25° x 0.25° grid. The region above declination 83° is 
left to a dedicated search for WIMP annihilation at the center 
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FIG. 10: Sky map of significances (— log 10 p) obtained in the full- 
sky search excluding trial factors (top), sky map of + v T 90% con- 
fidence level flux upper limits for an E~ 2 energy spectrum (10 -11 
TeV cm" 2 s _1 ) over the energy range 1.9 TeV to 2.5 PeV (middle), 
and the distribution of maximum significances for 1000 randomized 
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FIG. 11: Limits on an E~ 2 muon neutrino flux for sources in table 
|ni|from this work (triangles), limits from MACRO I14I1 . and Super- 
K fiUl . E~ 2 sensitivity for this work and the IceCube 22-string 
analysis, and predicted sensitivity for ANTARES 1 37] and IceCube. 
Our + v T limits are divided by 2 for comparison with limits on 
only Vy,. 



sky maps, with the obtained significance p ■ 
torn). 



7.4x10-" dotted (bot- 



' N by combining N sources of similar strength, with less 
improvement if one source is much stronger than average. We 
include five of eight sources and source candidates observed 
by Milagro with significance above 5a before considering trial 
factors, including four regions near Cygnus and one near the 
Equator. We add a hot spot near 8 = 1°, a = 19h [38], which 
may be associated with a large neutrino flux if confirmed as 
a source ITTIl . We exclude the three regions with pulsar-wind 
nebula counterparts, C3, C4, and the Crab Nebula, which are 
considered weaker candidates for significant hadron acceler- 
ation llllll . We adapt a method developed by HiRes [|39ll to 
perform our maximum likelihood search simultaneously for 
all six source locations, resulting in the slightly modified like- 
lihood function 



of the Earth l36ll . For each grid point, we similarly compute a 
flux limit and significance (Fig. [T0T >. We find a maximum pre- 
trial significance of p = 7.4x 1CT 4 at 6 = 54°, a = 11. 4h. We 
account for the trial factor associated with the all sky search 
by comparing the maximum pre-trial significance to the distri- 
bution of maximum pre-trial significances obtained from sky 
maps randomized in right ascension. We find 95% of sky 
maps randomized in right ascension have a maximum signifi- 
cance of at least p = 1 Ax 10 -4 (Fig. |T0b - Sensitivity and flux 
limits are summarized in Fig. QT| 

In the Northern Sky, the galactic TeV gamma ray sources 
observed by Milagro [2] are promising candidates for obser- 
vation with neutrino telescopes ifTH |12H. We improve our 
ability to detect a weak signal from this class of objects by 



i=l \ 3=1 / 

where S{ is the signal probability density of the i th event eval- 
uated for the j th source. Significance is again computed by 
comparing the obtained test statistic value to the distribution 
obtained from data randomized in right ascension. We ob- 
serve a small excess with a chance probability of 20%. The 
90% confidence level upper limit obtained on the mean flux 
per source is 9.7 x 10 -12 TeV cm" 2 s _1 . 

Finally, we search for groups of neutrino sources and ex- 
tended regions of neutrino emission by scanning for correla- 
tions of events at all angular distances up to 8°. We perform 
the search over a range of energy thresholds, using the num- 
ber of modules hit as an energy parameter. For each threshold 
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in angular distance and number of modules hit, we count the 
number of event pairs in the data and compare with the dis- 
tribution of pairs from data randomized in right ascension to 
compute significance. The highest obtained significance is p = 
0.1 with a threshold of 146 modules hit and 2.8° angular sep- 
aration, where we observe two event pairs. The probability of 
observing this maximum significance by chance is 99%. Since 
four separate analyses are performed on the data, the proba- 
bility of obtaining at least one significant result is increased. 
The most significant result obtained has a chance probability 
of 20%, and the binomial probability of obtaining this chance 
fluctuation in at least one of the four analyses is 59% and not 
significant. 

V. CONCLUSIONS 

We have analyzed 3.8 years of AMANDA-II data and found 
no evidence of high energy neutrino point sources. We place 
the most stringent limits to date on astrophysical point source 
fluxes. IceCube [40] is a next-generation neutrino telescope at 
the South Pole scheduled for completion in 201 1 with eighty 
60-module strings instrumenting ~1 km 3 of ice. Analysis of 
data recorded during 2007-2008 with the first 22 strings has 
improved the AMANDA-II sensitivity by a factor of 2. Cur- 
rently 59 strings are operating, and with continued construc- 



tion IceCube will achieve an angular resolution of better than 
one degree and an order of magnitude improvement over the 
AMANDA-II sensitivity within a few years. 
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